Efficient generation of random multipartite entangled states using time optimal 

unitary operations 
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We review the generation of random pure states using a protocol of repeated two qubit gates. We 
study the dependence of the convergence to states with Haar multipartite entanglement distribution. 
We investigate the optimal generation of such states in terms of the physical (real) time needed to 
apply the protocol, instead of the gate complexity point of view used in other works. This physical 
time can be obtained, for a given Hamiltonian, within the theoretical framework offered by the 
quantum brachistochrone formalism. Using an anisotropic Heisenberg Hamiltonian as an example, 
we find that different optimal quantum gates arise according to the optimality point of view used 
in each case. We also study how the convergence to random entangled states depends on different 
entanglement measures. 

PACS numbers: 03.67.-a; 03.67.Bg; 03.67.Ac 



CO 
CN 



43 ■ 
-t— > 



CN 
> 
O 

o 

CN 

in 
o 
oo 
o 



x 



I. INTRODUCTION 



One of the most fundamental concepts in the quan- 
tum description of Nature is that of entanglement [3, Q • 
Entanglement constitutes a physical resource that lies at 
the heart of important information processes @, 0, Q 
such as quantum tcleportation, superdense coding, and 
quantum computation. It has recently been proved that 
the generic entanglement, defined as the entanglement 
of random states, can be produced in a polynomial time 
using random two-qubit gates Q . 

Random states and random unitary operators are the 
quantum analogs of random numbers and are two basic 
concepts in quantum information and quantum commu- 
nication tasks. Random unitary operators are involved 
in the superdense coding of arbitrary states Q while the 
classical capacity of a noisy quantum channel is saturated 
by random quantum states [8|, just to mention two very 
significant applications. Both of them are well defined by 
the Haar measure which remains invariant under unitary 
transformation. Since the generation of random states 
is exponentially hard we study the production of states 
with the same distribution of entanglement, a task that 
can be performed with fewer physical resources Q. 

Although a large amount of different multipartite en- 
tanglement measures has recently been proposed, a con- 
siderable amount of research has particularly been de- 
voted to the study of multiqubit entanglement measures 
defined as the sum of bipartite entanglement measures 
over all (or an appropriate family of) the pos sible bi- 
partitions of the full system [H El, [ll, [H, [fl. There 
exist two popular entanglement measures for multiqubit 
pure states, one based on the von Neumann entropy of 
marginal density matrices and the other one based upon 
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the linear entropy of those matrices. It has recently been 
shown that the von Neumann entropy based measure is 
able to grasp more features of highly entangled states 
than the linear entropy based measure [l5| . It is then 
expected that when the efficient generation of random 
states is studied in the light of these entanglement mea- 
sures, some differences will arise because, as a conse- 
quence of the "concentration of measure" phenomenon, 
these states are almost maximally entangled (l6j . 

Emerson et al. [ItJ introduced a protocol for gener- 
ating pseudo-random unitary operators. A circuit of re- 
peated two-qubit gates was considered acting on a sepa- 
rable state of N qubits. The unitary operation is given 
by the combination of two independent single qubit rota- 
tions, chosen according to the invariant Haar measure at 
each time step, and a fixed two-qubit gate. The conver- 
gence to the Haar measure was shown to be polynomial 
with the number of qubits. Some experimental evidences 
have reinforced these results. A Markovian description of 
certain two-qubit gates has also been used to analytically 
prove that the convergence is reached in a polynomial 
time with the size of the system Q. Recently, a numer- 
ical effort has allowed to identify the optimal two-qubit 
gate improving this analytical bound [18j |. 

In these previous works convergence was studied in 
terms of gate complexity, i.e. the number of gates needed 
to reach the Haar distribution. In contraposition to the 
gate complexity, a new complexity concept for quan- 
tum al gori thms has been proposed: the time complex- 
ity (ill |20|. understood as the physical time needed to 
perform such algorithm. The minimization of this time 
is as important, from the experimental point of view, 
as the gate complexity. A considerable amount of work 
has recently been devoted to the time-optimal quantum 
computation problem, with emphasis in the quantum 
Brachistochrone formalism [2l|, [2^, [2^, [24[ ■ Making use 
of the analogy with the brachistochrone problem from 
classic mechanics, Carlini et. al [2l| introduced a varia- 
tional approach to obtain the optimal Hamiltonian and 
the optimal quantum evolution between initial and final 
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given states. A geometric approach to solve this problem, 
based on the symmetry properties of the quantum states 
space, was addressed in [22|. The role of the entangle- 
ment within the quantum brachistochrone formalism was 
studied in [H|. 

A more general result has recently been formulated in 
terms of the variational principle to find the time-optimal 
duration of a unitary transformation [2~4| . This formu- 
lation is independent of the input state and because of 
that more general than the previously described one. The 
time-optimal way to obtain a two-qubit universal quan- 
tum gate was previously discussed using the Cartan de- 
composition scheme for unitary transformation and un- 
der the constraint that one-qubit gates can be performed 
arbitrarily fast 0, [26J. To study the brachistochrone 
problem in relation with unitary operations, the quan- 
tum states space is replaced with the space of unitary 
operators. This formalism allows to consider the con- 
straint imposed by the finite amount of energy available 
in a physical experiment, as well as any other constraints 
imposed by experimental requirements or theoretical con- 
ditions. Then, the problem of finding the optimal param- 
eters for the Hamiltonian, is reduced to the resolution of 
a set of ordinary differential equations [24| . 

The aim of this contribution is to study the possible 
differences between both formalisms i.e. the gate and 
time complexities. We also study the dependence of the 
convergence rate with different entanglement measures. 
The paper is organized as follows. In Sec. II we review 
the protocol to generate random pure states and discuss 
the dependence of the convergence time, in terms of the 
number of gates to be applied, with the entanglement 
measure. In Sec. Ill we study the physical time for the 
convergence. Finally, Sec. IV is devoted to summarize 
and discuss our results. 



II. GENERATION OF RANDOM BIPARTITE 
ENTANGLEMENT 

The efficient quantum circuit generating random quan- 
tum states of N qubits is based on the iterative applica- 
tion of a two-qubit quantum gate Uij acting on qubits 
i and j, arbitrarily drawn from the N-qubit system, at 
each time step. The quantum gate Uij is composed by 
the product of two single qubit rotation gates Vi and Vj 
uniformly drawn from the Haar measure on U(2), and a 
fixed two qubit gate Wij 

U^ = ViVjW itj , (1) 
where Wij can be decomposed as 

Wij = (i>i®U2) exp [—% ^2 £>7c ® crfe] (ui®u 2 ), (2) 

k— x,y,z 

with the fixed rotations and u\,2 acting only on one 
of the two qubits, and are the Pauli matrices [27l. [28|. 



In our random entangling protocol the role of the fixed 
local rotations and should not be confused with 
the random local rotations Vij , which change at each 
step of the protocol. As we are just interested in the en- 
tanglement generation properties of the two-qubit gates 
we only need to consider the action of the non local part 
of the decomposition because its entangling power is 
the same when averaged over a large number of realiza- 
tions [28|, [29[ . The symmetries of such non-local action 
enables us to consider just a reduced range for the val- 
ues of its parameters (A& 6 [0,7r/4], k = x,y,z). The 
qubits i and j upon which the gate Uij is applied can be 
chosen in several different ways, and each of them corre- 
sponds to different geometries of the system: local and 
non local [30l ]. In the non local case qubits i and j are 
chosen randomly, the gate can act on two arbitrarily sep- 
arated qubits. In the local case the gate can only act on 
two neighboring qubits. In this scheme we study both, 
periodic and open boundary conditions. 

The results are qualitatively the same in the local and 
non local case, so we only show those of the non local 
couplings as representative of the typical behavior. The 
main difference is that a larger amount of two-qubit en- 
tangling gates are needed to converge to the entangle- 
ment of random states for the local geometry. 

Our goal is to reproduce the entanglement of typical 
random states, to such an end we averaged a large enough 
number of realizations (typically 10 3 ) in order to have 
small statistical fluctuations. 

The genuine multipartite entanglement E of a iV-qubit 
state can be expressed as 

1 [N/2] 

E = t— — - V E( m \ (3) 
[N/2] 2^ 

E (m) = JJm E ^(0- (4) 

bipart i=\ 

Here, E(i) stands for the entanglement associated with 
one, single bi-partition of the A-qubits system. The 
quantity gives the average entanglement between 

subsets of m qubits and the remaining N — m qubits 
constituting the system. The average is performed over 
the N^^ art nonequivalent ways to do such bi-partitions, 
which are given by 

K: P art = Q if N/2, (5) 

Kilt = 1(n /2 ) Xrn = N/2. (6) 

Different represent different entanglement proper- 

ties of the state, this is why all these entanglement mea- 
sures must be computed to capture all the entanglement 
properties of the state. The global multiqubit entangle- 
ment is given by the average of the [N/2] different £K m ) 
for any state 



3 




-13 B=B°f3°& D □ D -Ep 







E(« o 




1 E ° 






A A E<3 ' A " 












B □ A 















5 10 15 20 25 30 



Agates 



FIG. 1: Typical numerical simulation using the random cir- 
cuit in a system of 6-qubit. The entanglement average of the 
Haar measure, represented by dashed lines is reached in N ga tes 
steps. Top linear entropy as a measure of entanglement, bot- 
tom von Neumann entropy. We plot the global entanglement 
(squares) and the entanglement for the most balanced (trian- 
gles) and most unbalanced (circles) bipartitions of the system. 
Inset: Decay rates for the different entanglement measures, 
the rate is the same for the linear entropy and little differ- 
ences are observed for the von Neumann entropy for different 
bipartitions of the system. 



the average of all the single-qubit linear entropies [321 ] . 
This measure was later generalized by Scott to the case 
in which all possible bi-partitions of the system where 
considered [33j - 

We study the convergence according to different 
measures of multiqubit entanglement based upon bi- 
partitions. We characterized the global entanglement 
with the average of the bipartite entanglement measures 
associated with the 2^ — 1 bipartitions of the N-qubit 
system. 

We compare the evolution towards the convergence of 
the global entanglement with the evolution of the entan- 
glement of the most balanced bipartition E^ N / 2 ^ and 
the entanglement of the most unbalanced one fiW. We 
introduce the auxiliary normalized quantity 



AE 



Erq 



(E) 



(7) 



which decays exponentially with the number of itera- 
tions. AE will make easy the comparison between dif- 
ferent bipartitions. The saturation value Enaar, is the 
mean value of the entanglement of the Haar distribution 
given in [H, [34| and (E) is the averaged entanglement 
over system realizations. We choose the initial state to 
be the separable state |00 • • - 0) without loosing general- 
ity. The optimal time is defined as the number of gates 
required to reach AE = 0.01. We fixed the two-qubit 
gate A = (7r/4,0,0) which is known to be an optimal 
gate in the non-local couplings case [l|| . The results are 
qualitatively the same for any other choice of the gate. 

Fig. 1 shows that the convergence rates do not de- 
pend on the dimension of the chosen bipartition if the 
entanglement is quantified with the linear entropy How- 
ever, there exist small differences between the conver- 
gence rates if we use the von Neumann entropy as the 
measure of entanglement. We have a little faster conver- 
gence for the most unbalanced bipartition. This behavior 
is more visible for higher dimensional states, the larger 
the number of qubits, the larger the difference between 
the convergence rates. These results imply that when 
working with the linear entropy, it is enough to consider 
the convergence of the Meyer- Wallach entanglement mea- 
sure. If the von Neumann entanglement measure is used, 
one should consider the convergence of the entanglement 
of the most balanced bipartitions E^ N ^ 2 ^, because it is 
the one with a longer convergence time. 



We use two types of entanglement measures, El and 
E v n , respectively based on two different measures for the 
mixedness of the marginal density matrices pi associated 
with the bi-partitions: 

(i) the linear entropy Sl = 2 ™-i — Tr[p-]) and, 

(ii) the von Neumann entropy S v jy = — —Tr[pilogpi\. 

If one uses the linear entropy Sl, E^ turns out to 
be the well known Meyer- Wallach multipartite entangle- 
ment measure [3~H that Brennen showed to coincide with 



III. GATE COMPLEXITY VS. TIME 
COMPLEXITY 

In this section we are interested in the relation be- 
tween the number of gates in the circuit in order to reach 
the Haar distribution and the physical time to perform 
this operation. As an example we choose the Hcisenberg 
Hamiltonian: 

# = -E J ^k 2 + E s<v - ( 8 ) 
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where Jj are anisotropic couplings (j = x,y,z), B a (t) 
(a = 1,2) is an external magnetic field in z direction, 
and <rj = <7j (g> /, cr| — I <E) (Jj, with aj, as before, the 
Pauli matrices. The local magnetic terms appearing in 
([8]) are needed to perform the two-qubit gate operation 
and are not related to the random single qubit rotations 
Vij introduced in Sec. II. 

The optimal unitary evolution operator for the Hcisen- 
berg Hamiltonian (jSJ) was obtained in [24[ . Using the re- 
sults of the quantum brachistochrone formalism the op- 
timal time for some particular gates was calculated. 

We focus our efforts in the optimal entangler gate U$: 

(cos</> sin0\ 
— sin cos </>/ 

with the angle (f> G [0, n]. This gate applied to an initial 
separable state produces a </> dependent entangled state. 
When this gate is applied to the separable state |00), 
produces a maximally entangled state. The optimal time 
duration to implement the entangler gate assuming the 
finite energy condition is given by tot^ = 7ry / x(l — x/2), 
where x — ft/tr and cu is a constant given by the con- 
straint. Since the only effect of the external parameter 
cu is to rescale the time, in the following we take cu = 1 
as this choice will not modify our results. We study the 
convergence rate to the Haar distribution of this gate in 
terms of the physical time. This physical time is obtained 
as the product of the total number of gates N gate s times 
the optimal time {t p h ys = N gate s x t$). As in this 
case the final result is not affected by the choice of the 
measure of entanglement, we determine the convergence 
using El as a measure for the mixedness of the marginal 
density matrices. As before we consider the converge is 
reached for AE = 0.01 

Fig. 2 shows that a gate as simple as U<p, which 
depends just on the parameter <j), is enough to reveal 
the differences between the gate-complexity and time- 
complexity concepts. The curve N ga tes can be viewed 
as the time that would be necessary to reach the con- 
vergence if all the gates took the same time t<j, = 1 to 
perform each evolution. In such case, all gates with val- 
ues of (f> in the interval (7r/4, 37r/4) would contribute to a 
reasonably efficient algorithm. The optimal gate would 
be attained by <f> = 7r/2. But, as can be seen in the inset 
of Fig. 2, the optimal time needed to perform each 
iteration is far from being the same for all the values of 
<p. It increases with the value of cj>, been negligible for 
values of <j) near zero, point at which the entangling gate 
becomes the identity. If we combine both magnitudes, we 
obtain the total physical time t p h ys , which inherits the 
main properties of N gate s and t$. The global behavior 
resembles that of N gates , both of them diverge for the 
extremal values <f> = and (f> = n while attaining their 
lower and optimum values for central values of </>. The 
asymmetry of tphys compared to Ng ates comes from the 
behavior of t^. While N gat es is symmetric respect its 
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FIG. 2: Number of gates (dashed line) and physical time 
(solid line) to reach the convergence to the mean value of 
the Haar entanglement as a function of the gate parameter <j> 
for a system of 4 qubits. Inset: Optimal time to implement 
the entangler U$ gate as a function of <f>. As we have chosen 
units in which Planck's constant h is equal to one, all depicted 
quantities are dimensionless. 



optimum and central value <f> = 7r/2, the optimum gate 
according to t p h ys corresponds to a lower value of such 
parameter (0 ss tt/3). 

The case studied in this section is a good example 
of the role that the time complexity can play when de- 
signing an efficient quantum algorithm According to the 
gate complexity there exist a big range of values of (j) for 
which the random states generating algorithm is quite ef- 
ficient. If one introduces the time complexity argument 
this degeneracy is broken. The difference between the 
time needed by these efficient gates is not huge but it 
is enough to be taken into account, specially in a situa- 
tion where the algorithm must be run a large number of 
times. 



IV. SUMMARY AND DISCUSSION 

We have studied the generation of multipartite entan- 
gled states considering a protocol of a two-qubit fixed 
gate combined with two one-qubit random rotations. A 
comparison was made between the gate complexity and 
the time complexity revealing that they are two differ- 
ent optimality problems. In a real case both complexi- 
ties should be taken into account, and the optimal gate 
would be one with a reasonable good behavior according 
to both points of view. Nevertheless, when the algorithm 
must be run a large number of time, the optimization of 
time complexity is mandatory. In the example studied in 
Sec. Ill, the time complexity optimization allowed us to 
find an optimal gate between the whole family of gates 
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which arc almost equally efficient according to the gate 
complexity. 

The quantum brachistochronc formalism seems to be 
a promising approach for the treatment of the time com- 
plexity problem. Following the work of Carlini et al. [24[ 
it is possible to obtain the optimal way to realize a given 
quantum gate, introducing in this derivation any con- 
straint given by the experimental setup. 

We also focused in the possible dependence of the re- 
sulting entanglement on the different allowed bipartitions 
of the system. We found that when the linear entropy 
is considered, the convergence rates are independent of 
the bipartition scheme. However the convergence rates 
are different if we use the von Neumann entropy. In this 



case the most balanced bipartition should be considered 
to guarantee the convergence. These results arc indepen- 
dent of the possible, local or non local, geometries of the 
system. 
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